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Finite-size and surface effects in fine particle systems are investigated by Monte Carlo simulation of 
a model of a 7-Fe2C-3 (maghemite) single particle. Periodic boundary conditions have been used to 
simulate the bulk properties and the results compared with those for a spherical shaped particle with 
free boundaries to evidence the role played by the surface on the anomalous magnetic properties 
displayed by these systems at low temperatures. Several outcomes of the model are in qualitative 
agreement with the experimental findings. A reduction of the magnetic ordering temperature, 
spontaneous magnetization, and coercive field is observed as the particle size is decreased. Moreover, 
the hysteresis loops become elongated with high values of the differential susceptibility, resembling 
those from frustrated or disordered systems. These facts are consequence of the formation of a 
surface layer with higher degree of magnetic disorder than the core, which, for small sizes, dominates 
the magnetization processes of the particle. However, in contradiction with the assumptions of 
some authors, our model does not predict the freezing of the surface layer into a spin-glass-like 
state. The results indicate that magnetic disorder at the surface simply facilitates the thermal 
demagnetization of the particle at zero field, while the magnetization is increased at moderate 
fields, since surface disorder diminishes ferrimagnetic correlations within the particle. The change 
in shape of the hysteresis loops with the particle size demonstrates that the reversal mode is strongly 
influenced by the reduced atomic coordination and disorder at the surface. 

PACS Numbers: 05.10 Ln, 75.40 Cx, 75.40.Mg, 75.50 Gg, 75.50 Tf, 75.60 Ej 



I. INTRODUCTION 

The magnetic properties of fine particles are strongly 
influenced by finite-size and surface effects, their rele- 
vance increasing as the particle size decreases. Finite- 
size effects are due to the nanometric size of the par- 
ticles, while surface effects are related to the symme- 
try breaking of the crystal structure at the boundary of 
each particle. These effects are manifested in nanomet- 
ric particles through a wide variety of anomalous mag- 
netic properties with respect to those of bulk materials. 
The magnetic characterization of these systems has put 
forward the controversial issue of distinguishing between 
the contributions coming from finite-size and surface ef- 
fects to their peculiar magnetic properties. For instance, 
alternately explanations to the reduction of the satura- 
tion magnetization M s - a common experimental obser- 
vation in fine particle systems - has been given in the 
past. Early models postulated the existence of a so- 
called dead magnetic layer induced by the demagneti- 
zation of the surface spins, which causes a reduction in 
M s because of its paramagnetic response. In more re- 
cent works devoted to the study of different ferrimag- 
netic oxides - 7-Fe203, NiFe204, (CoFe204, CuFe20/t, 
in the form of nanometric particlesaB - a random spin 
canting at the surface, caused by competing antiferro- 
magnetic (AE) interactions, was observed by Mossbauer 
spectroscopy^!, polarizecB and inelasticDjieutron scatter- 
ing, and ferromagnetic (FM) resonanceQ. The origin of 
this non-collinear arrangement of the spins was [discussed 
by several authors sujpwcpking the surfaceoHiiEl or the 
finite-size explanationsliinLS, but up to the moment no 
clear conclusions have been established. 



All these ferrimagnetic fine particles share a singu- 
lar phenomenology at low temperatures. Among the 
static properties, experiments have shown that the hys- 
teresis kwjipsi— display high closure fields and do not 
saturatetSEJOeven at fields of the order of 50 T, which 
indicates that the anisotropy fields cannot be the only re- 
sponsible mechanism for the magnetization reversal. Low 
magnetization as compared to bulk, shifted loops after 
field cooling and irreversibilities between the field cool- 
ing and zero field,|Copljng processes even at high fields 
are also observed IjiMH. Moreover, the existence of ag- 
ing phenomenall3ll3 in the time-dependence of the mag- 
netization, indicates that there must be some kind of 
freezing leading to a complex hierarchy of energy levels. 
Whether these phenomena can be ascribed to intrinsic 
properties of the particle itself (spin-glass state of the sur- 
face wh ich creates an exchange field on the core of the 
particleliSli3) , or they are due to a collective behaviour 
induced by interpartiele interactionsE3~E3, has been the 
object of controversyEn 

Up to the moment there has been no model giving 
a clear-cut explanation of all the above mentioned phe- 
nomenology, but some works addressing part of the is- 
sues have been published in recent years. The first 
atomic-scale model of the magnetic behaviour of individ- 
ual ferrimaffnetic nanoparticles is due to Kodama and 
Berkowitzcil The authors presented results of calcula- 
tions of a micromagnetic model of maghemite particles 
which were based on an energy minimization procedure 
instead of the Monte Carlo (MC) method. They used 
Heisenberg spins with enhanced anisotropy at the sur- 
face with respect to the core and included vacancies and 
broken bonds at the surface, arguing that these are in- 
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deed necessary to obtain hysteresis loops with enhanced 
coercpity-^and high-field irreversibility. Later, Kachkachi 
et al.E3~E3 performed MC simulations of a maghemite 
particle described by a Heisenberg model, including ex- 
change and dipolar interactions, using surface exchange 
and anisotropy constants different to those of the bulk. 
Their study was mainly focused on the thermal variation 
of the surface (for them consisting of a shell of constant 
thickness) and core magnetization, concluding that sur- 
face anisotropy is responsible for the non-saturation of 
the magnetization at low temperatures. No attention 
was paid, however, to the magnetic properties under a 
magnetic field. 

Other computer simulations studying finite-size and 
surface effects on ferro- and antiferromagnetic Pqbic lat- 



tices have also been published. Trohidou et alJ£3'cfl per- 
formed MC simulations of AF small spherical clusters. 
By using an Ising model on a cubic latticeEa, they com- 
puted the thermal and magnetic field dependencies of 
the magnetization and structure factor, concluding that 
the particle behaved as a hollow magnetic shell. By 
means of a Heisenberg modeled with enhanced surface 
anisotropy, they studied the influence of different kinds 
of surface anisotropy on the magnetization reversal mech- 
anisms and on the temperature dependence of the switch- 
ing field. Dimitrov and WysinErEil studied the hysteresis 
phenomena of very small spherical and cubic FM fee clus- 
ters of Heisenberg spins by solving the Landau-Lifshitz 
equations. They observed an increase of the coercivity 
with decreasing cluster size and steps in the loops due to 
the reversal of surface spins at different fields. However 
they did not considered the finite temperature effects. 

In order to contribute to elucidate the above mentioned 
experimental controversies and to further develop the 
previously published numerical simulations, we present 
the results of a MC simulation of a single spherical par- 
ticle which aim at clarifying what is the specific role of 
the finite size and surface on the magnetic properties of 
the particle, disregarding the interparticle interactions 
effects. In particular, we will study the magnetic prop- 
erties under a magnetic field and at finite temperature, 
thus extending other simulation works. In choosing the 
model, we have tried to capture the main features of real 
particles with the minimum ingredients allowing to inter- 
pret the results without any other blurring effects. 

The rest of the article is organized as follows. In Sec. 
U we present the model of a maghemite particle upon 
which the MC simulations are based. In Sec. El the 
study of the basic equilibrium magnitudes - energy, spe- 
cific heat, and magnetization - in absence of magnetic 
field is presented, comparing results for different particle 



sizes with those for periodic boundaries. Sec. IV is de- 
voted to the study of magnetization processes under the 
presence of a magnetic field. The thermal dependence of 
hysteresis loops and coercive field are computed, and a 
detailed analysis of these quantities in terms of the sur- 
face and core contributions is performed. The effects of 
the introduction of different kinds of disorder on the mag- 



netic properties are presented in Sec. where we study 
both the deviation from ideal stoichiometry by random 
removal of magnetic ions on the whole particle, as well as 
the introduction of vacancies only at the surface of the 
particle. In Sec. VI, we end up with a discussion of the 
obtained results and a presentation of the conclusions. 



II. MODEL 

7-Fe203, maghemite, is one p£ the most commonly 
studied nanoparticle compoundscJ presenting the above 
mentioned phenomenology. Maghemite is a ferrimagnetic 
spinel in which the magnetic Fe 3+ ions with spin S = 5/2 
are disposed in two sublattices with different coordina- 
tion with the O 2- ions. Each unit cell (see Fig. [l]) has 8 
tetrahedric (T), 16 octahedric (O) sites, and one sixth of 
the O sites has randomly distributed vacancies to achieve 
neutrality charge. The T sublattice has larger coordina- 
tion than O, thus, while the spins in the T sublattice 
have Ntt = 4 nearest neighbours in T and Nto = 12 
in O, the spins in the O sublattice have Noo = 6 near- 
est neighbours in O and Not = 6 in T. In our model, 
the Fe 3+ magnetic ions are represented by Ising spins Sf 
distributed in two sublattices a — T, O of linear size N 
unit cells, thus the total number of spin sites is (247V 3 ). 
The choice of Ising spins allows to reproduce a case with 
strong uniaxial anisotropy, while keeping computational 
efforts within reasonable limits. The possible existence 
of a spin-glass state at the surface of the particle should 
be better checked by using an Ising model than one with 
continuous spins, since in the former frustration effects 
are enhanced^. Moreover, the Heisenberg version of the 
particle without disorder does not show irreversibility, 
the hysteresis loops, whereas the Ising version doesEill 
being easier to observe independently the effects of dis- 
order and finite size in the last case. 

The spins interact via antiferromagnetic (AF) ex- 
change interactions with the nearest neighbours on both 
sublattices and with an external magnetic field H, the 
corresponding Hamiltonian of the model being 

N a N al3 

n/k B = - J2 EE^ S " S H- 

a.fi=T.O i=l n=l 



a=T,0 i=l 



(1) 



where we have defined the field in temperature units as 

fiH 



h = 



(2) 



being S and /i the spin value and magnetic moment of 
the Fe 3+ ion, respectively. Hereafter, Si — ±1 and the 
maghemite values of the, nearest neighbour exchange con- 
stants will be consideredESEl: Jtt = -21 K, Jqo = -8.6 
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K, Jto = —28.1 K. Since the intersublattice interactions 
are stronger than those inside each sublattice, at low tem- 
peratures, there must be bulk ferrimagnetic order with 
spins in each sublattice ferromagnetically aligned and an- 
tiparallcl intrasublattice alignment. 

We have used periodic boundary conditions to simu- 
late the bulk properties and free boundaries for a spheri- 
cally shaped particle with D unit cells in diameter, when 
studying finite-size effects. In the latter case, two dif- 
ferent regions are distinguished in the particle: the sur- 
face formed by the outermost unit cells and an inter- 
nal core of diameter Dcore unit cells (see Fig. ||). The 
quantities measured after each MC step are the energy, 
specific heat, susceptibility and different magnetizations: 
sublattice magnetizations (Mo, Mr), surface and core 
magnetization [Msurfi Mc re)j an d total magnetization 
(MTotai)- Each of them have been normalized to the re- 
spective number of spins so that they can range from 1 
to -1. In particular, M^otai is 1 f° r ferromagnetic order, 
for a disordered system and 1/3 for ferrimagnetic order 
of the O and T sublattices. 

The size of the studied particles ranges from D = 3 
to 10 corresponding to real particle diameters from 25 to 
83 A (see Table |lj). In this table, we have also included the 
number of surface and core spins Ns U rf, Ncore, together 
with the normalized magnetization values of a ferrimag- 
netic configuration Mjj nc . Note that due to the finite size 
of the particles, the ratio of T and O spins produces Mjj nc 
values that, in general, do not coincide exactly with the 
1/3 value for perfect ferrimagnetic order in an infinite 
lattice. In order to make the measured magnetizations 
for different diameters comparable, we have normalized 
them to the corresponding Mjj nc values given in Table Q. 

III. EQUILIBRIUM PROPERTIES 

A. Energy and specific heat 

We start by studying the effect of free boundary condi- 
tions and finite-size effects on the equilibrium properties 
in zero magnetic field. The simulations have been per- 
formed using the standard Metropolis algorithm. Start- 
ing from a high enough temperature (T = 200 K) and an 
initially disordered state with spins randomly oriented, 
the system was cooled down at a constant temperature 
step ST = -2 K and, after discarding the first 1000 MC 
steps in order to allow the system to thermalize, the ther- 
mal averages of the thermodynamic quantities were com- 
puted at each temperature during a number of MC steps 
ranging from 10000 to 50000 depending on the system 
size. The starting configuration at each new temperature 
was the one obtained at the end of the averaging process 
at the previous temperature. Systems with periodic (PB) 
and free boundary (FB) conditions with spherical shape 
have been considered with sizes ranging from 3 to 14. 



In Fig. |3j, we compare the thermal dependence of the 
energy for spherical particles of different diameters D 
with the corresponding results for a system of size N = 14 
and PB (lowermost curve, left triangles). A second or- 
der transition from paramagnetic to ferrimagnetic order 
signaled by a sharp peak at T C (D) in the specific heat 
(see the Inset in Fig. ||) is clearly observed. Finite size 
effects on both the energy and the specific heat are very 
important even for D's as large as 14 in the FB case, 
while for PB conditions they are negligible already for 
N = 8. The energy difference between the disordered 
and ferrimagnetic phases as well as the critical tempera- 
ture T C (D) increases as D is increased. This last quantity 
is strongly size dependent and approaches the infinite size 
limit (T c (oo) = 126 ± 1 K as evaluated for the N = 14 
system with PB conditions) as D increases (see Fig. |^, 
in which the variation of the peak in the specific heat 
with 1/D has been plotted). T C {D) can be fitted to the 
scaling law 

r e M-r c (P) = (D_y 1/v 

T c (oo) \D J w 

as predicted by finite-size scaling theorjBS with Dq = 
1.86 ± 0.03 a microscopic length scale (in this case, it 
is roughly twice the cell parameter), and a critical ex- 
ponent v — 0.4£U± 0.03, which seems to indicate a mean 
field behaviouiE9. This result can be ascribed to the high 
coordination of the O and T sublattices. The fitted curve 
is drawn in Fig. ^| where deviations from scaling are ap- 
preciable for the smallest diameters for which corrections 
to the finite-size scaling of Eq. || may be importanta. 

B. Magnetization 

To study the effects of a free surface and of finite size on 
the magnetization of the particles, we compare in Fig. |E| 
the results for four particle diameters (D = 3, 4, 6, 8, open 
circles) with that corresponding to a N = 14 system with 
PB (representing the behaviour of the bulk). In this fig- 
ure, we have distinguished the surface (dashed lines) and 
core (dot-dashed lines) contributions to the total magne- 
tization (symbols). The results have been recorded dur- 
ing the same cooling procedure used to obtain the energy. 
The main feature observed is the reduction of the total 
magnetization Mxotai with respect to the PB case (con- 
tinuous line) due to the lower coordination of the spins 
at the surface, which hinders ferrimagnetic order at fi- 
nite temperatures. Fig. || clearly shows the roles played 
by the surface and the core in establishing the magnetic 
order. On one hand, independently of the particle size, 
the core (dot-dashed lines) tends to a perfect ferrimag- 
netic order at low T (marked by M = 1/3), progressively 
departing from the bulk behaviour as T approaches T c , 
this finite-size effect being more important as the par- 
ticle size decreases. However, the surface magnetization 
does not attain perfect ferrimagnetic order at T = even 
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for D — 8 due to the reduced coordination of the spins. 
For this reason, a rapid thermal demagnetization is ob- 
served which significantly departs Msurf from the bulk 
behaviour. 

It is worthwhile to note that for all the diameters stud- 
ied there is a temperature range in which this demagneti- 
zation process is linear, this range being wider as the par- 
ticle size decreases. In this linear regime, the particle de- 
magnetization becomes dominated by the surface effects, 
being the core and surface behaviours strongly correlated. 
Linear demagnetization is indicative of the effective 3D- 
2D dimensional reduction of the surface shell|-apd has 
previously been observed in thin files systemsEZrEa and 
in simulations of rough FM surfaces^. Mrotai is always 
strongly dominated by the surface contribution, progres- 
sively tending to the bulk behaviour as the particle size 
is increased. 



IV. HYSTERESIS LOOPS 

In Fig. ||, we show the hysteresis loops of particles with 
diameters D = 3, 6 for different temperatures. The loops 
have been computed by starting from a demagnetized 
state at h = and increasing the magnetic field in con- 
stant steps, 5h — 1 K, during which the magnetization 
was averaged over rj 3000 MC steps after thermaliza- 
tion. The results shown have been averaged for several 
independent runs starting with different random seeds. 

First of all, let us note that the saturation field and 
the high field susceptibility increase as the particle size 
is reduced, since this quantities are mainly associated to 
the progressive alignment of the surface spins towards the 
field direction. Thus, the loops of the smallest particles 
resemble those found in ferrimagnetic nap«Barticlesmi5cil 
and other bulk systems with disordeiE3c3, increasing 
their squaredness (associated to a uniform reversal mech- 
anism of M) with the size. In fact, by plotting sepa- 
ratedly the contributions of the core and the surface to 
the total magnetization (see Fig. [?], dashed lines), we see 
that the loop of the core is almost perfectly squared inde- 
pendently of temperature and particle size, indicating a 
uniform reversal of its magnetization with a well-defined 
ferrimagnetic moment. Instead, the loop of the surface 
reveals a progressive reversal of M, which is a typicaiica- 
ture associated to disordered or frustrated systemsEBElI 
Nonetheless, for a wide range of temperatures and par- 
ticle sizes, it is the reversal of the surface spins which 
triggers the reversal of the core. This is indicated by the 
fact that the coercive field of the core is slightly higher 
but very similar to the one of the surface. 

Since for all the studied particle sizes the h c (T) curves 
show a complex behaviour mainly related to the frus- 
tration of the antiferromagnetic intra and intersublatticc 
exchange interactions, we start by studying the case of a 
ferromagnet with no frustration. In Fig. the h c (T) 
dependence for a system with N = 8 and with the same 



lattice structure as maghemite but equal FM interactions 
Ja/3 = J and PB conditions is shown. The h c (T) de- 
pendence is now a monotonously decreasing curve with 
no inflection point, which at high enough temperatures 
(T/J > 1) can be fitted to a power law of the kind 

h c (T) = h c (0)[l-(T/T c ) 1 / a ] , (4) 

with a = 2.26±0.03; close but different to what would be 
obtained hy. a model of uniform reversal such as Stoner- 
Wohlfarthcd (a = 2). Even in this simple case, for which 
M reverses uniformly, the thermal variation of h c {T) can- 
not be only ascribed to the thermal activation of a con- 
stant magnetization vector over an energy barrier land- 
scape, since actually M is of course temperature depen- 
dent. Therefore, the reversal mechanism cannot be in- 
ferred from the a value obtained from a fit to Eq. || 
in any range of temperatures for which M significantly 
varies with T. 

The thermal dependence of h c for the maghemite par- 
ticles with AF interactions is shown in Fig. |j>. Both for 
the PB and spherical cases, the h c (T) curves are quali- 
tatively different from the FM case: they have opposite 
curvature and two regimes of thermal variation. 

Let us start by analyzing the PB case. At high T 
(T > 20 K), h c (T) can be fitted to the power law of Eq. 
| with a = 0.94 ± 0.02, h c (0) = 134 ± 2 K. Values of a 
close to 1 have been deduced in the past for some models 
of domain wall motioned At low T, a different regime is 
entered but tending to the same h c (Q) = 134.2 K. This 
change in behaviour is associated to the wandering of 
the system through metastable states with Mxotai — 0, 
which are induced by the frustration among AF inter- 
actions. Consequently, when lowering T, the hysteresis 
loops become step- like around h c as shown i n Fig j. b | (sim- 
ilar features are observed in related studiesEaE2lc3). The 
jumps at T = are located at h = 117, 134.2, 151.4 K, 
the values at which the magnetic field energy is enough 
to invert one O spin .having 0, 1, 2 O nearest neighbours 
inverted, respectivelyCJ. While the O sublattice does not 
reverse uniformly, the T sublattice instead reverses as a 
whole after the reversal of O, at ft, = 151.4 K. When T is 
increased from 0, the steps are rounded by the progres- 
sive population of states with greater degree of configu- 
rational disorder and less metastability, finally giving rise 
to the suppression of the steps for T around 12 K, when 
h c < 117 K and the high T regime of h c (T) is entered. 

The general h c (T) behaviour for spherical particles 
with FB strongly depends on the particle size. For 
D — 3,6 and T > 20 K, the h c decay is similar to that 
for PB, but, at any given T, being smaller than for PB 
and as the size of the particle is decreased. At these 
temperatures, h c (T) is dominated by the surface, which 
nucleates the reversal of the magnetization, as indicated 
by the proximity between the surface and core h c (see 
Fig. |7|). However, when lowering T below 20 K, Ms U rf 
and Mcore tend to be equal, the surface becomes less 
efficient as nucleation center for spin reversal, and h c be- 
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comes dominated by the core (h^ ur ^ < h^ ore for any 
particle size , see Fig. 0). This is the cause of the rapid 
increase of h c towards the PB values for the D — 6 curve 
(see Fig. |^b). For D = 3, instead, h c saturates when low- 
ering T due to the smaller ratio of core to surface spins, 
which actually hinders the prevalence of the core. 

Finally, it is worth noticing that, independently of the 
size of the particles with FB, the h c values are always 
smaller than that for PB, since the existence of spins 
with less coordination at the surface favours the forma- 
tion of reversed nuclei of spins acting as a seed for the 
reversal process, which is not the case of PB, where all 
equivalent spins have the same coordination. Therefore, 
the h c values for PB are only recovered at low T in the 
limit of large particle size, at difference with other exten- 
sive magnitudes such as the energy or the magnetization, 
for which we have checked that finite-size scaling is ac- 
complished. 



V. EFFECTS OF DISORDER 

In real particles, disorder and imperfections are present 
departing the system from perfect stoichiometry and dis- 
tort the position of the atoms on the .lattice, being their 
effect more important at the surfaceEj. There are sev- 
eral ways to implement this disorder on the model. The 
simplest way to simulate the deviation of the O and T 
sublattice atoms from ideal stoichiometry is by random 
removal of magnetic ions on the O/T sublattices. 



A. Disorder on the lattice 

Up to the moment, the existence of vacancies in the O 
sublattice in real maghemite structure has not been con- 
sidered. It is important to note that, in this system, intra 
and intersublattice magnetic interactions are antiferro- 
magnetic. Consequently, inclusion of vacancies in one of 
the sublattices may destabilize the FM parallel alignment 
of the other one, resulting in a system with a great de- 
gree of magnetic disorder. In particular, this effect will 
be much stronger when vacancies are introduced in the O 
sublattice, since Not is greater than Nto- To show the 
effect of these kind of disorder, we have simulated the 
hysteresis loops for different vacancy concentrations p v 
on the O sublattice at two cooling fields hpc — 20, 100 
K. As can be seen in Fig. [ToJ the introduction of a low 
concentration of vacancies (p v = 1/6 as in the real mate- 
rial) results in a reduction of the magnetization and in- 
creases the high field susceptibility without any substan- 
tial change in the general shape of the loops. However, 
if p v is increased beyond the actual value, the loops pro- 



B. Surface disorder 

In what follows, we will study the effects of the disor- 
der at the surface of the particle, considering a, p v = 1/6 
vacancy density on the O sublattice. Since the surface of 
the particles is not an ideal sphere, the outermost unit 
cells may have an increased number of vacancies on both 
sublattices with respect to those present in the core. Re- 
duced coordination at the surface may also change the 
number of links between the surface atoms. We will be 
denote by p sv the concentration of surface vacancies in 
the outermost primitive cells. 



1. Field Coolings 

The magnetic ordering of the system can be character- 
ized by studying the behaviour of the equilibrium magne- 
tization in a magnetic field. These curves have been ob- 
tained by the same cooling procedure used in the magne- 
tization simulations at zero field with ST = — 2 K in pres- 
ence of different cooling fields hpc- Several such curves 
are shown in Figs. [0], [l^, in which the surface (contin- 
uous lines) and the core (dashed lines) contributions to 
the total magnetization (open symbols) have been dis- 
tinguished. Let us first analize the case with no surface 
disorder (p sv = 0). The curves at different cooling fields 
do not collapse to the perfect ferrimagnetic order value at 
low T (i.e. M To tai = 1/4 for p sl) = l/6), reaching higher 
values of the magnetization the higher hpc, being this ef- 
fect greater as the particle size is reduced (compare Fig. 
|Tl| a and Fig. |l2]a). This is in contrast with the results 
for PB (not shown), for which the system reaches per- 
fect ferrimagnetic order at low T, even at fields higher 
than 100 K, evidencing that the main effect of the sur- 
face is the breaking of ferrimagnetic correlations within 
the particle. As a consequence, at a given temperature, 
the FM order induced by a magnetic field increases when 
decreasing D. 

By separatedly analyzing in detail the behaviour of 
the surface and core contributions to the total magne- 
tization, deeper understanding of finite-size effects can 
be gained. As in the case of ft, = 0, the total magne- 
tization for small particles is completely dominated by 
the surface contribution (continuous lines in Fig. 
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gressively closes, loosing squarec 
resembling those for a disorderec 
values of the high field susceptibilities and much lower 
coercivity. 



and progressively 
system, with high 



and this is the reason why the ferrimagnetic order is less 
perfect at these small sizes and the magnetic field can 
easily magnetize the system. However, the behaviour of 
the core of the smallest particles is still very similar to 
that of the case with PB, although its contribution to 
Mxotai is very small. At low fields, the surface is always 
in a more disordered state than the core: its magneti- 
zation lies below Mxotai at temperatures for which the 
thermal energy dominates the Zeeman energy of the field 
(see the continuous lines in Fig. [Tl|a, |l2|a). In this regime, 
the total magnetization closely follows that of the surface 
(see the curves in Figs. Ill and [l2| a for hpc — 20 K) 
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for the two simulated sizes. In contrast, a high field is 
able to magnetize the surface easier than the core due to 
the fact that the broken links at the surface worsen the 
fcrrimagnetic order, while the core spins align towards 
the field direction in a more coherent way. Only for the 
biggest particles the surface contribution departs from 
the Mxotai indicating the increasing contribution of the 
core (see the curves in Figs. |IT| a and fl2| a for hpc = 20 
K). Note also, that in this high Hfc regime, a maxi- 
mum appears which is due to the competition between 
the FM alignment induced by the field and the sponta- 
neous ferrimagnetic order (as the temperature is reduced 
the strength of the field is not enough as to reverse the 
spins into the field direction). 

The introduction of vacancies does not change the low 
field behaviour of the total magnetization, which is still 
dominated by the surface both for D = 3, 6, although 
the smallest particles are easily magnetized by the field. 
However, at high fields, Mxotai is lower than Mg ur f, the 
surface progressively decouples from M to tai with the in- 
troduction of vacancies in the surface, being this effect 
more remarkable for the biggest particle. With respect 
to the core, at difference with the non disordered case 
(Pv = Psv = 0), the low temperature plateau of Mcore 
tends to a higher value than that for perfect ferrimag- 
netic order, since the main effect of the disorder is to 
break ferrimagnetic correlations in the core; increasing 
the ferromagnetic order induced by the field. This is re- 
flected in a progressive departure of the high and low field 
Mcore curves with increasing disorder (see the dashed 
lines in the sequence b-d of Figs. |ll|, [l2]). The maxi- 
mum appearing at high /ipc is only slightly affected by 
disorder, shifting to lower temperatures and eventually 
disappearing for D = 3 and p sv = 0.5. 

2. Hysteresis loops 

Hysteresis loops with surface disorder are given in Fig. 
|l3| for two particle diameters. The introduction of sur- 
face vacancies facilitates the magnetization reversal by 
progressive rotation, producing a rounding of the hys- 
teresis loops when approaching h c , in the same way that 
occurs when particle size is reduced. The same fact ex- 
plains the increase of the high field susceptibility, since 
the vacancies act as nucleation centers of FM domains at 
the surface, which, from there on, extend the FM corre- 
lations to the inner shells of spins. Moreover, a consider- 
able decrease of h c is observed. All these facts yield to a 
progressive elongation of the loops, g»W»£ loop shapes re- 
sembling those of disordered systemg23 ! L2l. Fig. [l4|, where 
the surface and core contributions are shown separatedly, 
clearly evidences that the increase of FM correlations at 
the surface, facilitated by the vacancies, induce FM order 
in the core. That is to say, Mcore follows the evolution of 
Msurf at moderate fields above h c , in contrast with the 
case with no surface vacancies (see Fig. |7|) where the core 



keeps the ferrimagnetic order for the same field range. 

VI. DISCUSSION AND CONCLUSIONS 

We have presented a simple model of a maghemite 
nanoparticle with the minimal ingredients necessary to 
faithfully reproduce the magnetic structure of the real 
material. The model has proven successful in reproduc- 
ing several key features present in ferrimagnetic nanopar- 
ticle systems: (1) the reduction of T c , spontaneous mag- 
netization Mxotai, and coercive field h c , for small sizes, 
as D decreases; (2) the increase, with the reduction of the 
particle size and with the increase of surface disorder, of 
the differential susceptibility and the elongation of the 
hysteresis loops in resemblance with those of frustrated 
systems; and (3) the existence of a surface layer with 
higher magnetic disorder than the core. Let us further 
comment these points in deeper detail. 

First of all, we find that T C (D) follows conventional 
finite-size scaling, discarding any important surface ef- 
fect on this quantity. Similar finite-size effects have been 
found in fine particlesC-3 of MnFe2C>4, but with a sur- 
prising increase of T C (D) as D decreases. However, the 
spontaneous magnetization MxataU at any temperature, 
follows a quasi- linear behaviour with 1/D, see Fig. 
indicating that the reduction of MTotai is simply propor- 
tional to the ratio of surface to core spins, so it is mainly 
a surface effect. Similar experimental behaviour has been 
found in 7-Fe203C3 and the above mentioned systemo. 

The h c (T) thermal decay for the spherical particles is 
in qualitative agreement with the experimental results for 
maghemite particles of sizes 9-10 nm shown in Fig. 4 of 
Rcf. [I?], taking into account that in real samples there are 
additional contributions coming from the blocking pro- 
cess associated to the particle size distribution (h c drops 
to zero above the blocking temperature). In both cases, 
the curvature of the h c (T) curve is similar, suggesting a 
non-uniform reversal of the magnetization, a point that 
is also confirmed by the shape of the hysteresis loops 
around h c . However, our model for spherical particles 
gives reduced coercivities with respect to the bulk (rep- 
resented by the PB case). A fact that is ipj^antcast with 
the enhancement observed experimentallyli5l3E3, and in- 
dicating that finite-size effect cannot cause it. Increased 
anisotropy at the surface may be the responsible for it. 
In any case, the model qualitatively reproduces the h c 
reduction with D for small sizes [see Fig. 1 in Ref. OJ , 
which may be indeed a finite-size effect. 

The M(T) and M(h) dependencies obtained in our 
simulation lead to the conclusion that, in spherical par- 
ticles, there is a surface layer with much higher degree of 
magnetic disorder than the core, which is the Ising ver- 
sion of the random canting of surface spins [Occurring in 
several fine particle with spinel structureaEl. As_oppo- 
site with the suggestion given by some authorsBB that 
below a certain freezing temperature the surface layer en- 
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ters a spin-glass-like state, our model does not give any 
indication of this phenomenology at any of the studied 
sizes and temperatures. Furthermore, the surface layer, 
by partially breaking the fcrrimagnetic correlations, di- 
minishes the zero-field Mxotai but, at the same time, en- 
hances Mxotai at moderate fields. Although the surface is 
easily thermally demagnetized and easily magnetized by 
the field than the core, it does not behave as a dead layer, 
since, at any T, it is magnetically coupled to the core. 
All these facts put forward that the surface has higher 
magnetic response than the core, excluding a spin-glass 
freezing. Moreover, we do not observe irreversibilities 
between field and zero-field cooled magnetization curves, 
which is a key signature that in the scope of our model, 
neither finite-size or surface effects, nor the inclusion of 
surface vacancies are enough to account for the postu- 
lated spin-glass-like state. 

Finally, let us mention that our model does not repro- 
duce the experimentally observed shift of the hysteresis 
loops under field cooling adduced as a prove ofihcijcx- 
istence of the spin-glass-like state at the surfaceoEj til. 
Only when hp c smaller than irreversibility fields h irr 
are used in the numerical experiment, hysteresis loops 
that are apparently shifted are obtained, which in fact 
are minor loops. In any case, the absence of this phe- 
nomenology is in agreement with the non-observation of 
a spin-glass-like state at the surface, indicating that other 
ad hoc ingredients must be included in the model. For 
instance, enhanced surface anisotropy or exchange con- 
stants at the surface different than at-the bulk, as is the 
case in exchange coupled multilayersoeH. 

Current work is under progress to elucidate the possi- 
ble influence of these new ingredients and of interparticle 
interactions. 
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FIG. 1. Unit cell of maghemite. The magnetic Fe 3+ ions 
occupying the two sublattices, in different coordination with 
the 2 ~ ions (white colour), are coloured in black (T sub- 
lattice, tetrahedric coordination) and in grey (O sublattice, 
octahedric coordination). 
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TABLE I. Characteristic parameters of some of the spherical particles simulated: particle diameter D in units of the lattice 
constant a, diameter of the corresponding real D Rea i, number of total spins N To tal, number of spins at the surface and in the 
core Nsurf,Ncore, and magnetization of the noncompensated spins Mjj nc = (No — NT)/NTotai- The data are for particles 
with no vacancies in the O sublattice. 



D 


D Rcal (A) 


Nsurf 


Ncoro 


AiSurf 
M Unc 


iviCoro 
iV1 Unc 


M Unc 


3 


25 


330 (95%) 


17 


0.285 


0.412 


0.291 


4 


33 


724 (87%) 


111 


0.337 


0.369 


0.341 


5 


41 


1246 (78%) 


347 


0.355 


0.291 


0.341 


6 


50 


1860 (69%) 


841 


0.350 


0.332 


0.344 


8 


66 


3748 (58%) 


2731 


0.345 


0.330 


0.338 


10 


83 


6485 (48%) 


12617 


0.329 


0.337 


0.333 








40 80 120 160 200 
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FIG. 3. Thermal dependence of the energy for different diameters D = 3,6, 8, 14 (from the uppermost curve) and periodic 
boundary conditions N = 14 (lowermost curve). Inset: Thermal dependence of the specific heat for the same cases (the periodic 
boundary case is drawn with a dashed line). 
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1/D 

FIG. 4. Particle size dependence of the transition temperature T c from paramagnetic to ferrimagnetic phases for spherical 
particles with FB. The displayed values have been obtained from the maximum in the specific heat. The continuous line is a 
fit to Eq. 3. 
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FIG. 5. Thermal dependence of the magnetization M obtained by progressive cooling from high T at a constant rate, 
ST = —2 K, and starting from a random configuration of spins. The results for four particle diameters are shown: D — 3 (a), 
D — 4 (b), D = 6 (c), and D = 8 (d). The contributions of the surface (dashed line) and core spins (dot-dashed line) have 
been distinguished from the total magnetization (circles). The results for PB conditions, in a system of linear size iV = 14, 
have also been included for comparison (continuous line). 
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FIG. 6. Temperature dependence of the hysteresis loops for particles of diameter D = 3 (a), D — 6 (b). The temperatures 
starting from the outermost loop are T = 0, 5, 20, 40, 60, 80, 100 K, except for case D — 3 in which T = 0, 5 K curves are omitted 
for clarity. 
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aD Core 

FIG. 2. Schematic drawing of the spherical particles simu- 
lated in this study, showing the basic geometric parameters. 
The unit cells are indicated by the dashed grid, being the cell 
parameter a, and N the number of unit cells along each axis. 
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T(K) 

FIG. 8. (a) Temperature dependence of the coercive field h c for a system with the same structure as maghemite but 
ferromagnetic interactions (J a g = J) and PB conditions and N = 8; (b) Temperature dependence of the coercive field h c for 
the real AF values of the exchange constants for maghemite for the case of FB spherical particles of diameters D — 3 (circles), 
D = 6 (squares), and for a system of linear size N = 8 with PB conditions (diamonds). 
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FIG. 9. Detail of low temperature hysteresis loops for PB conditions around the coercive field h c . The corresponding 
temperatures are T = (circles), 5 K (squares), 10 K (diamonds), 20 K (triangles). 
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FIG. 10. Hysteresis loops for systems with vacancy concentrations p v = 0.0, 0.166, 0.4, 0.6 (from outer to innermost) on 
the O sublattice at T — 20 K. Particle diameters D = 3 (a), and D = 6 (b). Results have been averaged over 10 disorder 
realizations. 
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T (K) T (K) 

FIG. 11. Thermal dependence of M after cooling under a magnetic field for a spherical particle with D = 3, with vacancy 
densities on the surface of the O and T sublattices p av = (a), 0.1 (b), 0.2 (c), 0.5 (d), and p v = 0.166 on the O sublattice.The 
results for two cooling fields h,FC= 20, 100 K (lower and upper curves respectively in each pannel) are shown. The contributions 
of the surface (thick lines) and the core (dashed lines) to the total magnetization (circles) have been plotted separatedly. The 
magnetization has been normalized to Mb, the magnetization of a perfect ferrimagnetic configuration for a system of infinite 
size. 
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FIG. 13. Hysteresis loops for systems with vacancy densities on the surface of the O and T sublattices p sv ~ 0, 0.1, 0.2, 0.5, 
vacancy density p sv = 0.1666 on the O sublattice, and T = 20 K. Particle diameters D — 3 (a), D = 6 (b). Results have been 
averaged over 10 disorder realizations. 
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